MOSFiT: Modular open source fitter for transients

Much of the progress made in time-domain astronomy is accomplished by relating observational multi-wavelength time series data to models derived from our understanding of physical laws. This goal is typically accomplished by dividing the task in two: collecting data (observing), and constructing models to represent that data (theorizing). Owing to the natural tendency for specialization, a disconnect can develop between the best available theories and the best available data, potentially delaying advances in our understanding new classes of transients. We introduce MOSFiT: the Modular Open-Source Fitter for Transients, a Python-based package that downloads transient datasets from open online catalogs (e.g., the Open Supernova Catalog), generates Monte Carlo ensembles of semi-analytical light curve fits to those datasets and their associated Bayesian parameter posteriors, and optionally delivers the fitting results back to those same catalogs to make them available to the rest of the community. MOSFiT is designed to help bridge the gap between observations and theory in time-domain astronomy; in addition to making the application of existing models and creation of new models as simple as possible, MOSFiT yields statistically robust predictions for transient characteristics, with a standard output format that includes all the setup information necessary to reproduce a given result. As large-scale surveys such as LSST discover entirely new classes of transients, tools such as MOSFiT will be critical for enabling rapid comparison of models against data in statistically consistent, reproducible, and scientifically beneficial ways

Radio Emission from Ultra-Cool Dwarfs

The 2001 discovery of radio emission from ultra-cool dwarfs (UCDs), the very low-mass stars and brown dwarfs with spectral types of ~M7 and later, revealed that these objects can generate and dissipate powerful magnetic fields. Radio observations provide unparalleled insight into UCD magnetism: detections extend to brown dwarfs with temperatures <1000 K, where no other observational probes are effective. The data reveal that UCDs can generate strong (kG) fields, sometimes with a stable dipolar structure; that they can produce and retain nonthermal plasmas with electron acceleration extending to MeV energies; and that they can drive auroral current systems resulting in significant atmospheric energy deposition and powerful, coherent radio bursts. Still to be understood are the underlying dynamo processes, the precise means by which particles are accelerated around these objects, the observed diversity of magnetic phenomenologies, and how all of these factors change as the mass of the central object approaches that of Jupiter. The answers to these questions are doubly important because UCDs are both potential exoplanet hosts, as in the TRAPPIST-1 system, and analogues of extrasolar giant planets themselves.Comment: 19 pages; submitted chapter to the Handbook of Exoplanets, eds. Hans J. Deeg and Juan Antonio Belmonte (Springer-Verlag

An Unexpectedly Small Emission Region Size Inferred from Strong High-frequency Diffractive Scintillation in GRB 161219B

We present Karl G. Jansky Very Large Array radio observations of the long gamma-ray burst GRB 161219B (z = 0.147) spanning 1–37 GHz. The data exhibit unusual behavior, including sharp spectral peaks and minutes-timescale large-amplitude variability centered at 20 GHz and spanning the full frequency range. We attribute this behavior to scattering of the radio emission by the turbulent ionized Galactic interstellar medium (ISM), including both diffractive and refractive scintillation. However, the scintillation is much stronger than predicted by a model of the Galactic electron density distribution (NE2001); from the measured variability timescale and decorrelation bandwidth we infer a scattering measure of SM ≈ (8–70) × 10−4 kpc m−20/3 (up to 25 times larger than predicted in NE2001) and a scattering screen distance of d scr ≈ 0.2–3 kpc. We infer an emission region size of {\theta }_{s}\approx 0.9\mbox{--}4 μas (\approx (1\mbox{--}4)\times {10}^{16} cm) at ≈4 days, and find that prior to 8 days the source size is an order of magnitude smaller than model predictions for a uniformly illuminated disk or limb-brightened ring, indicating a slightly off-axis viewing angle or significant substructure in the emission region. Simultaneous multi-hour broadband radio observations of future GRB afterglows will allow us to characterize the scintillation more completely, and hence to probe the observer viewing angle, the evolution of the jet Lorentz factor, the structure of the afterglow emission regions, and ISM turbulence at high Galactic latitudes

Mapping Exoplanets

The varied surfaces and atmospheres of planets make them interesting places to live, explore, and study from afar. Unfortunately, the great distance to exoplanets makes it impossible to resolve their disk with current or near-term technology. It is still possible, however, to deduce spatial inhomogeneities in exoplanets provided that different regions are visible at different times---this can be due to rotation, orbital motion, and occultations by a star, planet, or moon. Astronomers have so far constructed maps of thermal emission and albedo for short period giant planets. These maps constrain atmospheric dynamics and cloud patterns in exotic atmospheres. In the future, exo-cartography could yield surface maps of terrestrial planets, hinting at the geophysical and geochemical processes that shape them.Comment: Updated chapter for Handbook of Exoplanets, eds. Deeg & Belmonte. 17 pages, including 6 figures and 4 pages of reference

Enabling New ALMA Science with Improved Support for Time-Domain Observations

While the Atacama Large Millimeter/submillimeter Array (ALMA) is a uniquely powerful telescope, its impact in certain fields of astrophysics has been limited by observatory policies rather than the telescope's innate technical capabilities. In particular, several observatory policies present challenges for observations of variable, mobile, and/or transient sources --- collectively referred to here as "time-domain" observations. In this whitepaper we identify some of these policies, describe the scientific applications they impair, and suggest changes that would increase ALMA's science impact in Cycle 6 and beyond. Parties interested in time-domain science with ALMA are encouraged to join the ALMA Time-domain Special Interest Group (ATSIG) by signing up for the ATSIG mailing list at https://groups.google.com/group/alma-td-sig

Enabling New ALMA Science with Improved Support for Time-Domain Observations

While the Atacama Large Millimeter/submillimeter Array (ALMA) is a uniquely powerful telescope, its impact in certain fields of astrophysics has been limited by observatory policies rather than the telescope's innate technical capabilities. In particular, several observatory policies present challenges for observations of variable, mobile, and/or transient sources --- collectively referred to here as "time-domain" observations. In this whitepaper we identify some of these policies, describe the scientific applications they impair, and suggest changes that would increase ALMA's science impact in Cycle 6 and beyond. Parties interested in time-domain science with ALMA are encouraged to join the ALMA Time-domain Special Interest Group (ATSIG) by signing up for the ATSIG mailing list at https://groups.google.com/group/alma-td-sig

Detection of cosmic structures using the bispectrum phase. II. First results from application to cosmic reionization using the Hydrogen Epoch of Reionization Array

Characterizing the epoch of reionization (EoR) at $z\gtrsim 6$ via the redshifted 21 cm line of neutral Hydrogen (HI) is critical to modern astrophysics and cosmology, and thus a key science goal of many current and planned low-frequency radio telescopes. The primary challenge to detecting this signal is the overwhelmingly bright foreground emission at these frequencies, placing stringent requirements on the knowledge of the instruments and inaccuracies in analyses. Results from these experiments have largely been limited not by thermal sensitivity but by systematics, particularly caused by the inability to calibrate the instrument to high accuracy. The interferometric bispectrum phase is immune to antenna-based calibration and errors therein, and presents an independent alternative to detect the EoR HI fluctuations while largely avoiding calibration systematics. Here, we provide a demonstration of this technique on a subset of data from the Hydrogen Epoch of Reionization Array (HERA) to place approximate constraints on the IGM brightness temperature. From this limited data, at $z=7.7$ we infer "$1\sigma$" upper limits on the IGM brightness temperature to be $\le 316$ "pseudo" mK at $\kappa_\parallel=0.33\,\textrm{"pseudo"}\,h$ Mpc$^{-1}$ (data-limited) and $\le 1000$ "pseudo" mK at $\kappa_\parallel=0.875\,\textrm{"pseudo"}\,h$ Mpc$^{-1}$ (noise-limited). The "pseudo" units denote only an approximate and not an exact correspondence to the actual distance scales and brightness temperatures. By propagating models in parallel to the data analysis, we confirm that the dynamic range required to separate the cosmic HI signal from the foregrounds is similar to that in standard approaches, and the power spectrum of the bispectrum phase is still data-limited (at $\gtrsim 10^6$ dynamic range) indicating scope for further improvement in sensitivity as the array build-out continues

Detection of cosmic structures using the bispectrum phase. II. First results from application to cosmic reionization using the Hydrogen Epoch of Reionization Array

Characterizing the epoch of reionization (EoR) at $z\gtrsim 6$ via the redshifted 21 cm line of neutral Hydrogen (HI) is critical to modern astrophysics and cosmology, and thus a key science goal of many current and planned low-frequency radio telescopes. The primary challenge to detecting this signal is the overwhelmingly bright foreground emission at these frequencies, placing stringent requirements on the knowledge of the instruments and inaccuracies in analyses. Results from these experiments have largely been limited not by thermal sensitivity but by systematics, particularly caused by the inability to calibrate the instrument to high accuracy. The interferometric bispectrum phase is immune to antenna-based calibration and errors therein, and presents an independent alternative to detect the EoR HI fluctuations while largely avoiding calibration systematics. Here, we provide a demonstration of this technique on a subset of data from the Hydrogen Epoch of Reionization Array (HERA) to place approximate constraints on the IGM brightness temperature. From this limited data, at $z=7.7$ we infer "$1\sigma$" upper limits on the IGM brightness temperature to be $\le 316$ "pseudo" mK at $\kappa_\parallel=0.33\,\textrm{"pseudo"}\,h$ Mpc$^{-1}$ (data-limited) and $\le 1000$ "pseudo" mK at $\kappa_\parallel=0.875\,\textrm{"pseudo"}\,h$ Mpc$^{-1}$ (noise-limited). The "pseudo" units denote only an approximate and not an exact correspondence to the actual distance scales and brightness temperatures. By propagating models in parallel to the data analysis, we confirm that the dynamic range required to separate the cosmic HI signal from the foregrounds is similar to that in standard approaches, and the power spectrum of the bispectrum phase is still data-limited (at $\gtrsim 10^6$ dynamic range) indicating scope for further improvement in sensitivity as the array build-out continues

Foreground modelling via Gaussian process regression: An application to HERA data

The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz